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Abstract 

Consider that a linear time-invariant (LTI) plant is given and that we wish to design a stabilizing 
controller for it. Admissible controllers are LTI and must comply with a pre-selected sparsity pattern. 
The sparsity pattern is assumed to be quadratically invariant (QI) with respect to the plant, which, from 
prior results, guarantees that there is a convex parametrization of all admissible stabilizing controllers 
provided that an initial admissible stable stabilizing controller is provided. This paper addresses the 
previously unsolved problem of determining necessary and sufficient conditions for the existence of 
an admissible stabilizing controller. The main idea is to cast the existence of such a controller as 
the feasibility of an exact model-matching problem with stability restrictions, which can be tackled 
using existing methods. Furthermore, we show that, when it exists, the solution of the model-matching 
problem can be used to compute an admissible stabilizing controller. This method also leads to a convex 
parametrization that may be viewed as an extension of Youla's classical approach so as to incorporate 
sparsity constraints. Applications of this parametrization on the design of norm-optimal controllers via 
convex methods are also explored. An illustrative example is provided, and a special case is discussed 
for which the exact model matching problem has a unique and easily computable solution. 
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I. Introduction 

In this paper, we deal with the problem of output-feedback stabilization for linear time- 
invariant (LTI) plants using sparsity-constrained LTI controllers^. The sparsity constraints are 
specified by a binary matrix with the same number of rows and columns as the controller. More 
specifically, entries of the controller must be zero whenever the corresponding element of the 
constraint matrix is zero, and are unrestricted otherwise. 

A. Previous Results 

The convex parametrization [2] proposed by Youla, which spans all LTI controllers that 
stabilize a prescribed LTI plant, popularized the so-called factorization approach [|3l to the 
analysis and synthesis of LTI feedback systems. The methods proposed in flU cast the search 
space in a ring, which provides additional insight and tools rooted on algebraic methods. 
However useful in expressing the design of norm-optimal controllers as convex programs, Youla's 
parametrization does not allow for sparsity constraints on the controller. The recent work in |Q, 
[0, 0, (El partially bridges this gap by identifying properties that the sparsity pattern of the 
plant and the one imposed on the controller must satisfy so that a convex parametrization of 
all stabilizing controllers may exist. These recently discovered methods spring from invariance 
principles that are valid in the presence of what the authors define as funnel causality, and their 
validity extended to the more general class flU of quadratically invariant sparsity patterns ifTOl . 
ifTTfl. The invariance condition in IfTOl . [[TT]| can be readily checked via an algebraic test, which, if 
true, assures that if there exists a stable stabilizing controller that satisfies the sparsity constraint 
then the set of all sparsity-constrained stabilizing controllers admits a convex parametrization 
based on a modification of the one in [[T2||. Subsequent work lfT3l has provided another convex 
parametrization that is guaranteed to exist under quadratic invariance, provided that a stabilizing 
controller that satisfies the sparsity constraint is given, and unlike prior work is not required 
to be itself stable. It has also been shown recently |fl4l that quadratic invariance of the set of 
controllers is necessary for the existence of the convex parametrization proposed in iflOll . [fTTTl . 

'For an interpretation of sparsity constraints in terms of the interconnection structure of distributed controller, see (T] Section 
III B]. 
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B. Contributions of this paper: 

The main results of this paper are motivated by the following problem. 

Problem 1.1. Consider that an LTI plant and a commensurate quadratically invariant sparsity 
constraint are given. Is the plant stabilizable by an LTI controller that satisfies the sparsity 
constraint? If one exists then compute it and give a convex parametrization of all stabilizing 
sparsity-constrained controllers. 

For a given plant, in this paper we establish necessary and sufficient conditions for the existence 
of a stabilizing LTI controller, subject to pre-specified quadratically invariant sparsity constraints. 
If one exists then our analysis also provides a method to construct a stabilizing controller that 
respects the sparsity constraints. Since all existing convex parametrizations presuppose prior 
knowledge of a stabilizing sparsity constrained controller ifTOl . our results bridge an important 
gap in the design process. 

In our solution method, the necessary and sufficient conditions mentioned in Problem [TTJ 
are cast as the existence of a certain doubly coprime factorization [fT51 . lfT6l of the plant that 
has additional constraints on the factors. We show that determining when such a factorization 
exists, and if so computing one, is equivalent to solving an exact model-matching problem 
with stability restrictions ifTTll . We also give a convex Youla-like parametrization of the set of all 
sparsity constrained stabilizing controllers by imposing additional constraints on the Q-parameter 
that require that it satisfies a certain homogeneous system of linear equations over the field of 
transfer functions. Unlike prior parametrizations that require an initial stable stabilizing controller 
that satisfies the sparsity constraint, our Youla-like parametrization does not require an initial 
controller and it is valid even when the plant is non-strongly stabilizable. 

C. Paper organization: 

Including the introduction, this paper has six sections. Section [II] states definitions and prelim- 
inary results used throughout the paper, while Section [III] reviews the notation and state of the 
art on design of sparsity constrained controllers. The main results of this paper are in Section ITVl 
where we formulate the necessary and sufficient conditions for stabilizability as the existence 
of solutions to an exact model matching problem IfTTll . We also propose methods to compute a 
sparsity-constrained stabilizing controller, when one exists, along with a numerical example. In 
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Fig. 1. Standard unity feedback configuration 



addition, we present an associated convex parametrization of all stabilizing sparsity-constrained 
controllers that is obtained by imposing subspace constraints on Youla's parameter. These results 
are specialized in Section [V] to plants that admit a structured doubly coprime factorization that 
we denominate Input/Output Decoupled. We show that this special factorization may simplify 
the application of our results and provide additional insights. The paper ends with conclusions 
in Section [VB 

Comparison with prior publications by the authors: Some of the results presented here have 
been published in preliminary form in lfT8l and |fT9l . In particular, parts of Sections HTT1 and HVl 
have been discussed with less detail in lfT8l . The discussion in [18J assumes block partitioning 
of the matrices, while, in this paper, partitioning is assumed only in Section |Vj In contrast with 
[fT8l . where we provide two simple examples, in Section ITV-BI of this paper we provide a more 
involved example on how to construct a sparsity-constrained controller. An abridged version of 
Section |V] was discussed in [19J in which there was a technical flaw. More specifically, we later 
found that Lemma 3.1 of |fT9l is incorrect and a correct and detailed discussion is provided in 
Section |V] and Appendix-I of this paper. This paper also establishes a strong connection between 
the approaches of Sections [IV] and |V] 

II. Preliminaries 

We focus on the standard feedback configuration of Fig. [fl where G is an LTI plant and K is 
an LTI controller that are finite dimensional and operate in either continuous or discrete-time. 
Here, v\ and v 2 are the input disturbance and sensor noise, respectively. In addition, u and y 
are the control and measurable output vectors, respectively. 
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We adopt the following notation: 



ITM /\ J 


OCL Ul d.11 ICdl — IctllUlldl Uu-llSlCl lUllCllUlla. 


R(A) nX9 


Spt of 7i x n matrices with pn trips in 1R( Ai 

iJV, l V 7 1 / L s\ {_J 1 1 1 CI LI IvvJ ¥V 1 LAI Willi IvJ All IL^i. 1 s\ I • 


1 rlvl 


Transfer function matrix, or, equivalently, |^J n 


n 


Stability region for TFMs. 


A(A) 


Subset of R(A) whose entries have poles in Vt. 


A(A) nx <? 


Set of n x q stable TFMs. 


A(A) 


Set of stable TFMs, or, equivalently, {J n q A(\) nxq 


B 


The set {0,1}. 



We also adopt the following assumptions and conventions: 



m Dimension of y. 

p Dimension of u. 

G The plant is a TFM with strictly proper entries. 

K The LTI controller is an element of R(A) pxm . 
H(G, K) The TFM from [u[ vl f to [y T u T f. 

g) Kronecker product. 

~q {l,2,---,g} 



The indeterminate A is either s for continuous-time or z for discrete-time systems, respectively. 
The A argument of a TFM is often omitted when its presence is clear from the context. If the 
transfer matrix H(G, K) is stable we say that K is a stabilizing controller of G, or equivalently 
that K stabilizes G. If a stabilizing controller of G exists, we say that G is stabilizable. 



A. Coprime and Doubly Coprime Factorizations for LTI Systems 

A right coprime factorization (RCF) of G over Vt is a fractional representation of the form 
G = NM~\ with N G A mxp and M G A pxp , and for which there exist X G A pxm and 
Y G A pxp satisfying YM + XN = I ([3, Ch. 4, Corollary 17]). Analogously, a left coprime 
factorization (LCF) of G (over Q) is defined by G = M~ l N, with N G A mxp and M G A mxm , 
satisfying MY + NX = I for X G A pxm and Y G A mxm . Due to the natural interpretation 
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of the coprime factorizations as fractional representations, the invertible M and M factors are 
sometimes called the "denominator" TFMs of the coprime factorization. 

Definition 11.1. A collection of eight stable TFMs (M,N, M, N, X,Y, X,Y) is called a 
doubly coprime factorization (DCF) of G over Vt if M and M are invertible, yield the following 
factorizations: 

G = M~ X N = NM' 1 
and satisfy the following equality (Bezout's identity): 



Y 


X 




M 


-X 


-N 


M 




N 


Y 



<-m+p- 



(1) 



To avoid excessive terminology, we refer to doubly coprime factorizations over Vt simply as 
doubly coprime factorizations (DCFs) (31 Ch.4, Remark pp. 79]. 

Theorem II.2. (Youla) [3 Ch.5, Theorem 1] Let (M, N, M, N, X, Y, X, Y) be a DCF of G. 

Any stabilizing controller can be written as: 



K Q = XqYq 1 = Yq 1 Xq 
for some Q in A pxm , where Xq, Xq, Yq and Yq are defined as: 

X Q = f X + QM 
X Q = f X + MQ 
Y Q = f Y-QN 



Y 



def 



Q 



Y-NQ 



(2) 

(3) 
(4) 
(5) 
(6) 



It also holds that Kq stabilizes G for any Q in A pxm . 

Remark IL3. The following identity shows that (M, N, M, N, Xq, Yq, Xq, Yq) in Theorem \II.2\ 
is also a DCF of G: 

(Y-QN) (X + QM) M -{X + MQ) 
—N M N (Y-NQ) 



L m+pi 



Q e 



(7) 
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III. Feedback Control subject to Sparsity Constraints 

The precise formulation of the sparsity constrained stabilization problem is achieved by 
imposing a certain pre-selected sparsity pattern on the set of admissible stabilizing controllers. 

A. Specifications of Sparsity Constraints on LTI Controllers 

For the boolean algebra, the operations (+, •) are defined as usual: 0+0 = 0-1 = 1-0 = 0-0 = 
and 1 + = + 1 = 1 + 1 = 1-1 = 1. By a binary matrix we mean a matrix whose entries 
belong to the set B. With the usual extension of notation, B 9X ' stands for the set of all binary 
matrices with q rows and I columns. The addition and multiplication of binary matrices is carried 
out in the usual way, keeping in mind that the binary operations (+, •) follow the boolean 
algebra. Binary matrices are marked with a "bin" superscript, in order to distinguish them from 
transfer function matrices over R(A). Furthermore, for binary matrices of the same dimension, 
the notation A hm < B hm means that < bij holds entry wise for all i and j. 

A binary matrix may be associated with a TFM of the same dimension, whereby each entry 
of the binary matrix corresponds to an entry of the TFM. The following definitions introduce 
operators that will be used to establish a correspondence between binary matrices and the sparsity 
pattern of G or sparsity constraints imposed on K. 

Definition III.l. (Pattern operator) Given A in R(X) qxl , we define <p(A) e M qxl as follows: 

/ as d *f I if =0 . . 

p(A)ij — < i,jel,qxl,l (8) 

1 otherwise 

Definition III.2. (Sparse operator) Conversely, for any binary matrix A bm in M qxl , we define 
the following linear subspace: 

§(A bin ) = f {A g R(\) qxl \ p(A) < A bin | (9) 

Definition IIL3. Given K bm in W xm , the sparsity constraint S is defined as KWJ : 

S= f §(K hln ), (10) 

Hence, S is the subspace of all controllers K in M(A) pxm for which K tj = whenever = 0. 

We assume that S and p act on G and C7 bin on an analogous way as above, leading to the 
following definitions. 
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Definition III.4. The following is the sparsity pattern of G (G hin G M mx P): 

G hin = f p(G) (11) 

Remark III.5. From matrix multiplication (in the boolean algebra), we conclude that the 
following holds: 

p(K G) < p(K) p(G) (12) 

B. Quadratic Invariance 

Definition III.6. 470] Definition 2] The sparsity constraint S is called quadratically invariant 
(QI) under the plant G if 

KGK G S, K eS. (13) 

Remark III.7. The following conditions are equivalent to rfTil) / fiOl/ : 
. p(KGK) < K bin , K eS 

Definition III.8. Define the feedback transformation h G : M(A) pxm M(A) pxm of G with K, 
as follows: 

h G (K)= f K(I + GK)~\ KeR(\) pxm (14) 

Remark III.9. The feedback transformation hc(-) is invertible, and its inverse is given by 

h G \K) = f K(I - GK)' 1 (15) 

Proof: Note that h G {-) is well-posed because K is proper and G is strictly proper. The rest 
of the proof follows by direct algebraic computations and is omitted for brevity. ■ 
The following result is used throughout the paper. 

Theorem 111.10. / TiQl Theorem 14] Given a sparsity constraint S, the following equivalence 
holds: 

S is QI under G <^=^ h G (S) = S (16) 
where we adopt the following abuse of notation: 

h G (S) = f {h G {K)\K G S] 
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An alternative algebraic proof of Theorem (TT6l ) is given in [20, Theorem 9]. 

Remark III.ll. The set S is QI under the given plant G if and only if S is QI under —G. This 
implies, via ([76]) above, that S is QI under G if and only if Iiq 1 (S) = S, where Hq 1 (S) d = 
{h G \K)\KeS}. 

IV. Main Result 

Given a QI sparsity constraint S, in Theorem IIV.2I we develop necessary and sufficient 
conditions for the existence of a stabilizing controller in S. These conditions are formulated 
in terms of the existence of a doubly coprime factorization of the plant in which the factors 
satisfy additional constraints. Such a factorization (when it exists) is equivalent to solving an exact 
model matching problem with stability [17] restrictions, which has been previously investigated 
in the control literature. The following preparatory result will be used throughout this Section. 

Proposition IV.l. Let (M, N, M, N, X, Y, X, Y) be a given DCF of G. The following identities 
hold: 

MXq = Kq (I + GKq) -1 , XqM = Kq(I + GKq) 1 , Q e A pxm (17) 

Proof: We proceed to verifying that MXq = Kq(I+ GKq) 1 is true, while the proof that 
XqM = Kq(I + GKq) 1 holds is omitted because it is analogous. From Kq = Yq 1 Xq and 
G = NM~\ we get that K Q (l + GKq)' 1 = (I + Y^XqNM'^Y^Xq, where we used 
the fact that Kq (J + GKq) 1 = (/ + KqG) 1 Kq. Finally, using Bezout's identity we find that 
(I + Y^XqNM- 1 )' 1 = (I + Y Q \I -YqM)M- 1 )" 1 = MY Q , which by direct substitution in 
(/ + Y^XqNM- 1 ) ^Y Q l X Q concludes the proof. ■ 
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The following Theorem is a main result of this paper. 



Theorem IV.2. Let (M, N, M, N, X, Y, X, Y) be a DCF ofG and SbeaQI sparsity constraint. 

• Sufficiency: If Q in A pxm is such that at least one of the following inequalities hold^: 

p(X Q M) < K bin (18a) 
p(MX Q ) < K hin (18b) 

then Kq is a stabilizing controller in S. 

• Necessity: If there is a stabilizing controller in S then there exists some Q in A pxm for 
which both inequalities in rfTffl) hold and, in addition, the controller can be written as Kq. 



Proof: Necessity: Suppose that there exists a stabilizing controller in S, then, as a conse- 
quence of Youla's Theorem 111.21 such a controller can be written as Kq for some Q in A pxm . 
We now use the fact that Kq is in S to prove that both inequalities in ([TBI hold. According to 
Proposition IIV.1I we get from (fT7T) that 

X Q M = K Q (I + GKqY 1 (19) 

We apply the p operator © on both sides of equation (fT9l) and using Definition IIII.8I we find 
that p(XqM) = p{h G (K Q )). Since S is QI and K Q is in S, it follows from (OH) that h G (K Q ) 
belongs to S and p(h G (K Q )) < K bin , which leads to p(X Q M) < K bin . Similarly, we employ 
(fTTT) to get that p(MXq) = p(h G {K Q )) in order to finally obtain that p(MX Q ) < K bin . 

Sufficiency: Take each side of ([191 as an argument for h G l {-) in order to get via Definition IIII. 8 1 
that h G l (X Q M) = h^ihc^Ko)) and equivalently that K Q = h G 1 (X Q M). In addition, it follows 
from Remark IIII.9I and Remark IIII. 1 II that /i G; 1 (iS) = S, which, from the assumption that 
p(XqM) < K hin , implies that Kq = h G 1 (X Q M) is in S . The fact that K Q is stabilizing 
follows from Youla's Theorem 111.21 The sufficiency with respect to p(MX Q ) < K bm follows 
by a similar line of proof and so is omitted for brevity. ■ 



2 In fact, it also follows from the statement of the Theorem that either both inequalities hold or none. 



DRAFT 



11 

A. Controller Synthesis as An Exact Model-Matching Problem with Stability Restrictions 
Henceforth, given a matrix V with n rows and q columns, we adopt the following notation: 

vec(V) v = vec(V) gives % + (j_i )n ) = V id 
diag(V) A = diag(V) is diagonal and An = Vi 

In this section, we will outline a method (based on Theorem IIV.2I above) for the computation 
of a stabilizing controller subject to a pre-selected QI sparsity constraint S (whenever such a 
controller exists). Given a DCF of G, which can be computed using the standard state-space 
techniques in |fT5ll . lfT6ll . our goal is to obtain Q in A pxm such that (fT8l) is satisfied. 

Our approach is based on the realization that (fl"8T) can be cast as the feasibility of an exact 
model-matching problem ifTTl with respect to Q in A pxm . This correspondence is stated precisely 
in the following Theorem, while Section IIV-A1I provides more details and references on the 
computation and tests for the existence of solutions of the exact model matching problem. 



Theorem IV.3. Consider that a DCF of G (M, N, M, N, X,Y, X ,Y) is given and that a QI 
sparsity constraint S is pre-selected via a choice of K bm in B pxm . The existence of a stabilizing 
controller in S is equivalent to the existence of Q in A pxm for which at least one of the following 
equivalent equalities holds: 

$ (M T ® M) vec(Q) + $ vec (XM) = (20a) 
$(M T ® M)vec(Q) + $ vec (MX) = (20b) 

where $ is defined as: 

<$>= f I- diag(X bin ) (21) 
In addition, if there is Q in A pxm that satisfies d20l ) then Kq is a stabilizing controller in S. 



Proof: The proof follows by establishing an equivalence between (l20l) and (TT8T) . We start 
by rewriting (TT8T) as follows: 

p((X + MQ)M) = p(MQM + XM) < K hin (22) 
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p(M(X + QM)) = p(MQM + MX) < K' 



bin 



(23) 



The vectorization of (|22l) -(|23l leads to: 




1) Computational considerations: Problems of the type (|20l) are of particular importance 
in linear control theory and were formulated and proposed for the first time by Wolovich 
( ifTTlO . who also coined the term exact model-matching in the early 1970's. Under the additional 
constraint that Q lies in A pxm , the problem is referred to as exact model-matching with stability 
restrictions (see EQ). Reliable and efficient state-space algorithms for solving (|2Q|) are available 
in E2l . which also describes a method to ascertain when a solution exists and consequently, from 
Theorem IIV.31 decide when a stabilizing controller in S exists. Given a stabilizing controller 
in S one can use the results in lfT3l . [|23l to obtain a convex parametrization of all stabilizing 
controllers in S. Also, since the resulting convex parametrization is affine in Q, one can use the 
tractable methods proposed in ifTOll to design norm-optimal controllers for both the disturbance 
attenuation and the mixed-sensitivity I-L2 problems. 

B. A Numerical Example 

Consider the following choices for the plant G and the QI sparsity constraint S to be imposed 
on the controller as specified via K bm : 



A+4 



1 

A-2 



G 



A-1 







A+5 



1 

A-3 



bin 







1 



1 



1 1 



We use the state-space formulas from [fT5l . [fl~6ll to obtain the following DCF of G: 
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M 



A-2 
A+6 










A-l A-3 
A+7 A+7 

A-3 
A+8 







N 



A-2 



1 



(A+4)(A+6) A+6 

2(A+1) 1 

(A+5)(A+7) A+7 

A-3 1 

(A+8)(A+5) A+8 



X 



A-2 1 A-3 

A+6 A+7 A+8 

1 A-l 1 

A+6 A+7 A+8 



Y 



A 5 +42A 4 +617A 3 +4144A 2 +12968A+15256 A 4 +777A 3 +9557A 2 +27300A- 13060 
A 5 +30A 4 +355A 3 +2070A 2 +5944A+6720 

2A+14 



N 



A 2 +14A+48 
A-l A-3 



A 5 +30A 4 +355A 3 +2070A 2 +5944A+6720 
A 2 +40A-108.002 
A 2 +14A+48 



(A+4)(A+9) (A+10)(A+11) 

A+9 
A-l A-2 



M 



A-l 
A+9 







(A-2)(A-3) 
(A+10)(A+11) 



(A+5)(A+9) (A+10)(A+11) 

The remaining factors X and Y of the DCF are not needed here. We now proceed to finding a 
solution for (T20l) . which, according to Theorem IIV.31 leads to the conclusion that a stabilizing 
controller in S exists. In addition, we will use the aforementioned solution to compute a 
stabilizing controller. 

Since there are two zeros in the sparsity pattern imposed by K bm , the system of equations in 
(|2Q|) has two (nontrivial) equations that are satisfied by the following element of A pxm : 



Q 



1 o 1 

1 ^±§ 
U 1 A+7 

-1 



The resulting stabilizing central controller K = Yq Xq is given by 



K 



A+17 
A+7 

(A+5.87)(A-0.4525) 
' 04 (A+4)(A+5)(A+6)(A+8) 





(A+42.5389)(A-2.5389) 
(A+6) (A+8) 



AT7 

A+6 u A+8 



which is in S. 
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C. A Youla-like Parametrization of All Sparse, Stabilizing Controllers 

In this subsection, we present an alternative statement to Theorem IIV.3I that clarifies the 
differences between it and Youla's classical parametrization. 



Corollary IV.4. Let S be a given QI sparsity constraint and (M, N, M, N, X, Y, X, Y) a DCF 
of G. Assume that there is a stabilizing controller in S and let Qq in A pxm be selected to satisfy 
d20l ). Any stabilizing controller in S can be written as Kq, where Q is obtained as: 

Q = Qo + Qs (26) 

for some Q$ in the (convex set) specified by the following inclusions: 

vec(Q s ) e Null($(M T ®M)),Q 5 e A pxm (27) 

where $ is the matrix defined in ([27]). 



Proof: The proof follows directly from Theorem IIV.3I ■ 
Corollary IIV.4I unveils the fact that once any suitable Q is found then the set of all stabilizing 
controllers in S can be generated from the affine subspace specified by (|26l)-(|27T). Notice that in 
Youla's classical approach the parameter Q is only required to be in A pxm , while the additional 
constraints in (l26l)- (|27T) guarantee that the resulting controller will be in S. 

D. Numerical Considerations 

For an introduction to linear subspaces for TFMs and vector bases of such subspaces we 
refer to [24J. In addition, the authors of [1231 describe a systematic, state-space algorithm to 
determine a basis of the null space of $(M T cg> M). Note that the main result in [|25l enables 
the computation of a basis having only stable poles, by performing a column compression of 
the normal rank of $(M T eg) M) by post-multiplication with a unimodular matrix. Furthermore, 
this basis is also minimal, in the sense that the basis-matrix, obtained by juxtaposing the basis 
columns, has no Smith zeros. Hence, this may be used for the parametrization of all stable 
vec(Q) in Null($(M T g> M)). 

For a numerical example illustrating Corollary IIV.41 from Subsection IIV-CI above, we refer to 
Subsections IV-C and IV-E in ifTSl . 
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E. Norm-Optimal Control Design 

We now indicate how Theorem IIV.2I can be used in conjunction with results from |[T0l to 
design norm-optimal controllers. In particular, given a quadratically invariant sparsity constraint 
S one may be interested in solving the following optimization problem: 

min \\H{G,K)\\ (28) 

K stabilizing and K(^S 

where || • || is a suitably defined operatorial norm. 
Using Theorem IIV.2I we can rewrite (l28l) as follows: 

min_ \\H{G,K Q )\\ (29) 
(MX + MQM) G S 

Q e A pxm 

where we used the fact that the inequalities in (fT8~l) are equivalent to XM + MQM G S and 
MX + MQM G S. Notice that Theorem |IV3] guarantees that (EU) is feasible if and only if ([29]) 
is feasible. 

We proceed by noticing that the closed loop TFM for a given controller Kq can be written 

as: 

H(G,K Q ) 



YM ■ 


-YN 




N 


Q 












+ 




M 


N 


(30a) 


XM I 


-XN 




M 










I - NX 


-NY 




N 














+ 




Q 


M 


N 


(30b) 


MX 


MY 




M 











where we used the formulae available in [3, pp.110]. Hence, we can use (1291) to rewrite (1281) as 
follows: 



min_ ||Ti + T 2 QT 3 || (31) 
(MX + MQM) G S 

Q G A pxm 



where T l5 T 2 and T 3 are obtained from (1301) . 

The analysis above implies that the sparsity constrained disturbance attenuation problem (as 
introduced in ifTOl (l)/pp. 276 ]), or the sparsity constrained mixed H 2 sensitivity problem (from 
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11261 pp. 139]) can be solved as a model-matching problem via the numerical technique in [fTOl 
Theorem 29]. 

V. Block-decoupling and streamlined solutions 

In this Section, we consider that the input and the output vectors of G are partitioned into 
blocks so that, under certain conditions, G can be factored in a special form that simplifies both 
the solution of the exact model matching problem of Theorem HV.3I and the parametrization in 
Corollary |IV.4[ Henceforward, we consider the following notation: 



number of partitions of y 



{y[i]}i=i partitions of the output vector 
nii dimension of ym 

r u number of partitions of u 

{ u [f\Yi=i partitions of the input vector 
Pi dimension of U[q 

The partitions are constructed in a way that the following holds: 



u 



T T 
V[X] ' ' ' V[r v ] 



T T 
U [l] ' ' ' U [ru 



1 T 



m 



V 



Similarly, we also consider the partitioning of G and K as: 



G 



G[\l] ■ ■ ■ G[i ru ] 



G[ ry l] ■ ■ ■ G [ryr . u] 



K[lT 



K 



(32) 



(33) 



[r u r y ] 



DRAFT 



17 



Assumption V.l. Throughout this Section, we assume that G and an associated partition of the 
input and output 021) are given. 

Remark V.2. Given factorizations of G and K as G = M _1 iV = NM" 1 and K = XY' 1 = 
Y~ 1 X, respectively, the partition in ( TJ21 will induce a unique block-partition structure on the 
factors N, M, N, M, X, Y, X and Y as well. 

A. Input/Output Decoupled Coprime Factorizations 

We start by defining input and output decoupled factorizations for G. 

Definition V.3. Let N and M be a factorization of G. The pair (N, M) is called output decoupled 
if M has the following block diagonal structure: 



M = diag({M M }2L i ; 



(34) 



where diag({M[ ii ]}^ 1 ) is defined as: 



diagtt/%]}^) d = f 



Af[n] 
M [22l 
















Mr, 



(35) 



Definition V.4. Let N and M be a factorization of G. The pair (N, M) is called input decoupled 
if M has the following block diagonal structure: 



M = diag({M N }£ 1 ) 



(36) 



Remark V.5. Notice that an output decoupled factorization can always be constructed by 
factoring each block row of G separately as follows: 



G[n] ■ ■ ■ G [iru] ] = M,J 



Nu 



i e 1,ty 



(37) 



An input decoupled factorization can also be constructed by factoring the block columns of G. 

Definition V.6. A DCF (M, N, M, N, X, Y, X, Y) of G is called input/output decoupled if the 
pairs (N, M) and (N, M) are input and output decoupled, respectively. 
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It is important to notice that the procedure outlined in Remark |V5l does not guarantee that the 
pairs (N, M) and (N, M) will be co-prime, much less doubly co-prime. In fact, G may not admit 
an input/output decoupled DCF. Sufficient conditions and algorithms to obtain an input/output 
decoupled DCF for G are provided in the Appendix I. 

There are two substantial benefits of working with an input/output DCF for G: The first is that 
the constraint on Q in Theorem IIV.2I reduces to Q E S f] A pxm , which leads to a parametrization 
of all stabilizing controllers that has a simpler characterization. The second advantage is that 
the exact model-matching problem of Theorem IIV.3I admits a unique solution that can be easily 
computed (see Section IV-CI) . 

B. Theorem \IV.2\ Revisited 

Here, we modify the definitions of Section UU] so that they account for the assumed 
input/output partition in (|32|) . More specifically, sparsity constraints will be imposed on 
entire block sub-matrices of K. The definitions in Section [III] can be recovered from the ones 
below for the case when the block sub-matrices have dimension one, i.e., provided that r y = m 
and r u = p. 

Definition V.7. Given K in R(A) pxm , we define p(K) G l r « xr * as follows : 

(v\ def J ° if K M = Q PiXmi • - ^i x , i ,o C x 

p( K )ij = \ t,j el,r u x l,r y (38) 

I 1 otherwise 

where PiXmj is a matrix with rows and rrij columns and whose entries are all zero. 

Definition V.8. Conversely, for any binary matrix K hm in W uXry , we define the following linear 
subspace: 

S(K hin ) = f \k e R(X) pxm \p(K) < K hin } (39) 

Definition V.9. Given K hm in M TuXry , the sparsity constraint S is defined as: 

S= f §(K Un ), (40) 

Hence, S is the subspace of all controllers K in R(A) pxm for which = 0p iXmj whenever 
K^ n = 0. In addition, we assume that § and p act on G and G bm as well as on the factors of 
any DCF of G in an analogous way. 
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Remark V.10. As a consequence of the definitions above, the following holds for any input/output 
decoupled DCF of G: 

p(M) = I ruXru , p(N) < G bm 

~ ~ {a) (41) 

p(M) = I ryXry , p(N) < G bm 

(b) 

where (a)-(b) follow from (O and the fact that G = M~ l N = NM^ 1 . 

The following Corollary is an immediate consequence of Theorem IIV.2I and the facts that 
p(M) = I ru xr u , p(M) = Ir y xr y , and M- 1 and M" 1 are well defined TFMs. 

Corollary V.ll. Let (M, N, M } N, X,Y, X,Y) be an input/output decoupled DCF of G and 
S be a Ql sparsity constraint. 

• Sufficiency: If Q in A pxm is such that at least one of the following inequalities hold^: 

p(X Q ) < K bin (42a) 
p(X Q ) < K bin (42b) 

then Kq is a stabilizing controller in S. 

• Necessity: If there is a stabilizing controller in S then there exists some Q in A pxm for 
which both inequalities in j[42\) hold and, in addition, the controller can be written as Kq. 



The following Corollary is the main result of this section. 

Corollary V.12. Let S be a given QI sparsity constraint and (M, N, M, N, X, Y, X, Y) an 
input/output decoupled DCF of G. Assume that there is a stabilizing controller in S and let Qq 
in A pxm be selected to satisfy ft42\) . Any stabilizing controller in S can be written as Kq, where 
Q is obtained as: 

Q = Qo + Qs, Q s eSn A pxm (43) 

Proof: From Corollary IV. 1 1 1 and Theorem II V. 3 1 it follows that since Q satisfies (|42l) then it 
will also satisfy (1201) . Hence, from Corollary IIV.4I any stabilizing controller in S can be written 
as Kq, with Q = Qo + Qs, where Qs satisfies (l27l) . The proof follows from noticing that since 

3 In fact, it also follows from the statement of the Theorem that either both inequalities hold or none. 
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M g) M is block diagonal and its inverse is a well defined TFM, Q satisfies (1271) if and only if 
vec{Q s ) G Null($) n A pxm holds, or equivalently Q 5 is in 5 n A pxm . ■ 

C. Theorem \IV.3\ revisited 

In this subsection, we show that if an input/output decoupled DCF of G exists then we can 
use Corollary IV. 1 1 1 to obtain a simplified version of Theorem II V. 3 [ A precise statement of this 
result is given in Corollary IV.15I 

Definition V.13. We define the binary matrix K± n belonging to the set M TuXry as follows: 

. ... def f 1 if K% n = 0, 

K m ) l3 = { 13 (44) 

otherwise . 

Definition V.14. Given K^ n G M ruXTy we introduce the linear subspace S± ofR(X) pxm as 

S ± = f | if e R{X) pxrn p{K) < (45) 

Corollary V.15. Let (M, N, M, N, X, Y , X, F) be an input/output decoupled DCF of G. Given 
a QI sparsity constraint S, G is stabilizable by a controller in S if and only if M~ 1 X$ ± is in 
A pxm , where X$ x results from the additive factorization X = Xs+X$ ± satisfying p(X$) < K hm 
and p(X s± ) < K h l n . 



Proof: We start by restating the first equation of (1421) for any Q as: 

p(X + MQ) = p(X s + MQ S + M(M~ 1 X S± + Q s± )) < K bin (46) 

where Q = Q s + Q s± and p(Q s ) < K hm and p(Q s± ) < K^ n - 

We now recall that according to Corollary IV. 1 1 1 G is stabilizable by a controller in S if and 
only if there is Q in A pxm so that (l46l) holds. However, given the fact that M _1 is block diagonal, 
(l46l) holds for some Q in A pxm if and only if M~ 1 X$ ± = —Qs ± has a solution where Qs ± is 
in A pxm . Since M~ l X$^ = —Qs ± nas a unique solution because M is invertible, we conclude 
that there exists Q s± is in A pxm satisfying (l46l) if and only if M~ 1 X S± is in A pxm . Notice that 
p(M~ 1 X s± ) < K± n always holds because M is block diagonal. ■ 
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VI. Conclusions 

We address the design of stabilizing controllers subject to a pre-selected quadratically invariant 
sparsity pattern. We show that the previously unsolved problem of determining stabilizability 
with sparsity constraints is equivalent to the solvability of an exact model-matching system of 
equations that is tractable via existing techniques, and we also outline a systematic method to 
compute an admissible controller. The proposed analysis also leads to a convex parametrization 
that is an extension of Youla's classical result so as to incorporate sparsity constraints on the set 
of stabilizing controllers. We indicate how this parametrization can be used to write sparsity - 
constrained norm-optimal control problems in convex form. 

Appendix I 

This Appendix has two parts. In the first part we give a sufficient condition that guarantees 
the existence of an output (input) decoupled left (right) coprime factorization for G, as in 
Definitions IV. 3 1 and IV 41 In the second part, we show that if G admits the aforementioned 
factorizations then there is a state-space method for computing its input/output decoupled DCF 
of Definition IV6l 



A. A Sufficient Condition for the Existence of an Output (Input) Decoupled Left (Right) Coprime 
Factorization 

We are given a plant G, partitioned as in (|33l) . As described in Remark IV 5 1 we perform a 
left coprime factorization for each of the r y block-rows of G (such left coprime factorizations 



always exist and can be computed using the classical state-space methods from 111511 . 111611 ) in 
order to get 

[G [tl] ■ ■ ■ G [iru] ] = Mr 1 



N[n] ■ ■ ■ N[ ir 



iel,r v . (47) 



where the poles of Mm can be placed anywhere in the stability domain f2. 

The following proposition gives a necessary and sufficient condition under which the row 
factorizations in (|47l) can be concatenated to produce a left decoupled coprime factorization for 
G. It should be noted that a left decoupled coprime factorization for G may exist that cannot 
be constructed from the row factorizations in (1471) . This indicates that the proposition is only a 
sufficient condition for the existence of a left decoupled coprime factorization for G. 
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Proposition VI.l. Let (M,N) be an output decoupled left factorization derived from the 
coprime factorizations rt47l) as follows: 



row 



G 



[ii] 



G 



[!»•«] 



M [n, 







Mr 1 , 



an J consider \& to &e f/ze following TFM: 

M[ii] ••• iV [n] • 



[ii] 



AT, 



l'u 



(48) 



<fe/ 



lr u ] 



••• M M N [ryl] 



The following holds: 

1) The output decoupled left factorization in ft48\) is coprime if and only if the following holds: 



rank(\l>(A)) = m, A G A 



G 



(49) 



where A G represents the set of unstable poles of G. 
2) The condition in ft49\) does not depend on the choice of the row coprime factorizations 

m. 



Proof: The proof follows as a consequence of standard results in linear systems theory, so 
we only provide a sketch of the ideas. We start by invoking a result used in [27] that the left 
factorization (M, N) is coprime if and only if ^ has no Smith zeros^l in C — . Hence, from 
the statement in 1), we are left to prove that any Smith zero of ^ in C — is a pole of G. 
The proof of 1) is concluded by noticing that if a given A in C — f2 is not a pole of G then, 
from the coprimeness of the row factorizations in (|47l) , M (A ) is invertible and hence full rank, 
leading to the conclusion that ^(Aq) must also be full row rank, and hence A is not a Smith 
zero of ^. 

It only remains to prove 2). The argument here follows from the fact that the set of all left 
coprime factorizations of any block-row of G is given by (PTTl) up to a premultiplication by a 
unimodular TFM [fJl Ch. 4, Theorem 43]. This in turn implies that ^(A) is unique, up to a 



A complex number Ao G C is called a Smith zero of ^ if $(Ao) is not full-row rank. 



DRAFT 



23 

premultiplication of a block-diagonal, unimodular TFM (having the same block partition as M), 
which does not alter the rank condition on \I/(A) at any unstable A G (C — Vt). ■ 
The corresponding test for the existence of input decoupled right coprime factorization of G 
(Definition IV. 41) is analogous and is therefore omitted for brevity. 



B. A State-Space Method to Compute an Input/Output Decoupled DCF 

We assume that G admits an output decoupled left coprime factorization G = M _1 N and 
an input decoupled right coprime factorization G = NM~ X . These factorizations must be pre- 
computed (using for instance the arguments described in Appendix I- A), so we consider the 
M, M, N, N factors fixed. Under these conditions, here we provide a computational algorithm 
to obtain the Input/Output Decoupled DCF of G (Definition [V6]>, i.e. a DCF 



Y 


X 




M 


-X 


-N 


M 




N 


Y 



(50) 



containing the fixed factors M, M, N, N. The fact that the coprime factorization (|50l) always 
exists is guaranteed by (31 Ch. 4, Theorem 60] but since we are not aware of any method to 
actually compute it, we will present one here. 

The following additional notation is needed: given any n-dimensional state-space representa- 
tion (A, B, C, D) of a LTI system, its input-output description is given by the transfer function 
matrix (TFM) which is the m x p matrix with real, rational functions entries 



G 



A 


B 


C 


D 



def 



D + C(AJ n — A)~ l B, 



(51) 



where A, B,C, D are nxn, nxp, mxn, mxp real matrices, respectively while n is also 
called the order of the realization. For elementary notions in linear systems theory, such as 
controllability, observability, detectability, we refer to [|28l . or any other standard text book in 
linear systems. 

We have started out with an output decoupled left coprime factorization (Definition IV. 3 1) of 
the plant G = M~ l N . The state-space representation of this left coprime factorization can be 
obtained according to Proposition IVI.3I A) from Appendix II, starting from a certain stabilizable 
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state-space realization of G (which we take without loss of generality to be in the Kalman 
Structural Decomposition, ||29l ) and which we consider fixed: 



G 



k 


k 


k 


k 


k 





^22 





^24 


B 2 








-k 
























c 2 





c. 


D 



(52) 



with the * denoting parts of the realizations that are of no importance in this proof. Continuing 
with Proposition IVI.3I A) from Appendix II, there also exists an invertible matrix U and a 
feedback matrix F (both fixed) such that (with F partitioned in accordance with (l52l) ) we get 



—N M 



u- 1 



* 


k 


k 




k 




k 





^4-22 ~~ FiC% 





A24 


— F2C4 


B 2 - F 2 D 


F 2 





k 


k 




k 


-k 


k 





—F4C2 





Au 


— F4C4 


-F A D 


F, 





-c 2 







-c 4 


-D 


I 



(53) 



A 



c n. 



with 



A22 — F2C2 A24 — F2C4 
— F4C2 A44 — F4C4 

Note that since (|52l) is stabilizable it follows that K{A^) C fi. After removing the unobservable 
part from (l53l) we get that 



(54) 







A22 ~ F2C2 


A 24 


— F2C4 


B 2 - F 2 D 


F 2 


N M 


= u- 1 


— F4C2 


A M 


— F4C4 


-F4D 


F 4 






-C2 




-c 4 


-D 


I 



(55) 



We have also started out with an input decoupled right coprime factorization (Definition IV.4I) 
G = NM -1 . According to Proposition IVI.3I B) in Appendix II, there exists a certain detectable 
state-space realization of G (which we take without loss of generality to be in the Kalman 
Structural Decomposition) and which we also consider fixed: 
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' A n 


A 12 




* 









A 22 







B 2 


G = 










-k 
















★ 










c 2 





* 


D 



(56) 



with the * denoting parts of the realization that are of no importance here. Any two realizations 
of G will always have the same the controllable and observable part, up to a similarity trans- 
formation - that is to say that if the controllable and stabilizable part of (|52l ) is (A 22 , B 2 ,C 2 , D) 
then the controllable and stabilizable part of (l56l) must be (Z~ 1 A 22 Z,Z~ 1 B 2 , C 2 Z,D), for 
some invertible, real matrix Z . We can apply this similarity transformation adequately on (|56l) . 
such that the the controllable and stabilizable part (A 22 , B 2 ,C 2 , D), appears identical on both 
realizations (|52|) and (l56l) . respectively. This simplifies future computations. 

According to Proposition IVI.3I B) from Appendix II, along with realization (l56l) . there also 
exists an invertible matrix V and a feedback matrix L (both fixed) such that (with L partitioned 
in accordance with (l56l) ) 







' A n - B X L X 


A 12 


— B\L 2 




•k 


Bx 






—B 2 Li 


A 22 


— B 2 L 2 




* 


B 2 


M 














•k 





N 















•k 









-In 




-L 2 




-k 


I 






—DL\ 


c 2 


- DL 2 


* 


-k 


D 



with 



V 



(57) 



A 



c n, 



(58) 



An — B\L\ A\ 2 — B\L 2 
—B 2 L\ A 22 — B 2 L 2 

Note that since (l56l) is detectable it follows that A(A U ) C fi. After removing the uncontrollable 
part from (1571) we get that 
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' A u - B\L\ 


A12 


— BiL 2 


Bi 


M 




—B2L1 


A22 


— B2L2 


B 2 


N 




-L x 




-L 2 


I 






-DU 


Co 


~DL 2 


D 



V 



(59) 



We have come now to the following stabilizable and detectable state-space realization of G, 
which we consider fixed: 





'A u 


A 12 


* 


B x 


G = 





A22 


A24 


B 2 














A M 










c 2 


c 4 


D 



(60) 



Since A(A U ) C ft we get that ([60]) is detectable and since A(^4 44 ) C ft we get that ((60l) is 
stabilizable, hence (|60l ) satisfies the hypothesis from Theorem I VI . 21 iii) from Appendix II. Starting 
from realization (|60l ) (which is fixed), (|68l) and (|69l ) yield a valid DCF of G for any stabilizing 
feedback matrices F + and L + (partitioned in accordance with (l60l) and satisfying Theorem IVI.2I 
ii) from Appendix II), and any invertible matrix T + satisfying the hypothesis of Theorem I VI.2| i). 
We will denote the factors of this particular DCF with (M+, A^+, M + , A^+, X + , Y+, X + , Y+) . 
After removing the unobservable part, the M + factor will be (the computation are similar with 
those for getting from (O to (1551) ) 



where 





^4-22 — F2C2 


A 2 4 




c 4 




M + = U + 1 


-F+C2 


A44 




c 4 






-C 2 




-c 4 




I 



(61) 



A 



A 



22 



F^Co Ao 



-FiC 2 



2 W2 /124 

A M 



F 2 C 4 
FfC\ 



c n. 



(62) 



and U + is a real, invertible matrix. We compute the factor 6 d = MM+ using the state-space 
realizations from (1531) and (loTI) respectively and we get 
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Q = U 



-1 



A22 — F2C2 


A24 


— F2C4 


F2C2 


F2C4 


F 2 ' 


—F4C2 


A44 


— F4C4 


F4C2 


F4C4 


F4 










A22 


A24 


F + 

r 2 













A44 


F + 
r 4 


-c 2 




-Ci 


c 2 


C4 


I 



After removing the unobservable part from (1631) we get that 





A22 ~ F2C2 


a 2 4 


— F2C4 


F 2 


- F + 


e = u- 1 


—F4C2 


A44 


— F4C4 


F 4 


- F + 
1 4 




-c 2 




-C4 


I 



u + 



and consequently 





A22 ~ F2C2 


A24 


- F 2 + C 4 


F 2 


- F + 

r 2 




-F+C 2 


A44 


— F^C4 


F4 


- F + 
r 4 




c 2 




C4 


I 



u, 



(63) 



(64) 



(65) 



which combined with (|54l) and (l62l) shows that 6 is unimodular. A similar line of reasoning can 
be used to prove that d = M + ~ l M is unimodular. 
Finally, compute 



0- 1 O 
O 6 



Y+ X+ 
-N+ M+ 



M+ -X+ 
N+ Y + 



O 

o ©- 1 



'm+p 



(66) 



which is still a DCF of G in its own right, due to the unimodularity of and 0. Plugging in 
the definitions of and into (l66l) yields 



Q-ix+ 
-N M 



M -X+Q- 1 
N Y+Q- 1 



(67) 



which is an input/output decoupled DCF of G and the algorithm ends. 
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Appendix II 

Theorem VI.2. [30 Theorem I] Let G be some proper m x p TFM. The class of all DCFs ([7]) 
of G over Q is given by 



M -X 
N Y 



Y X 
-N M 



T -i 



A-BL 


B 


F 






-L 


I 





T 




C-DL 


D 


I 






A-FC 


B 


- FD 


F 


L 




I 







-C 




-D 




I 



(68) 



(69) 



where A, B,C, D, F, L and T are real matrices accordingly dimensioned such that 
V W 



i)T 



O U 



has its diagonal p x p block V and m x m block U respectively, invertible, 



ii) F and L are feedback-matrices such that A(A — BL) U A(A — FC) C Q, 



Hi) G 



A 


B 


C 


D 



is a stabilizable and detectable realization. 



Corollary VI.3. Let G be an arbitrary m x p TFM and Q a domain in C. 

A) The class of all left coprime factorizations of G over Q, G — M~ X N, is given by 



N M 



U~ 



1 


A-FC 


B-FD -F 




C 


D I 



(70) 



where A, B,C, D, F and U are real matrices accordingly dimensioned such that 

i) U is any m x m invertible matrix, 

ii) F is any feedback matrix that allocates the observable modes of the (C, A) pair to Q, 



Hi) G 



A 


B 


C 


D 



is a stabilizable realization. 



B) The class of all right coprime factorizations of G over Q, G = NM is given by 



M 
N 



A 


- BL 


B 






-L 


I 


V 


C 


- DL 


D 





(71) 



where A, B,C, D, L and V are real matrices accordingly dimensioned such that 
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i) V is any p x p invertible matrix, 

ii) L is any feedback matrix that allocates the controllable modes of the (A, B) pair to Q, 



Hi) G 



A 


B 


C 


D 



is a detectable realization. 



The proof of Corollary IVI.3I follows on the lines of Theorem IVI.2I 
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